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Abstract 

We analyze the properties of deuteron-like structures in infinite, correlated 
nuclear matter, described by a realistic hamiltonian containing the Urbana f 14 
two-nucleon and the Urbana TNI many-body potentials. The distribution of 
neutron-proton pairs, carrying the deuteron quantum numbers, is obtained as 
a function of the total momentum by computing the overlap between the nu- 
clear matter in its ground state and the deuteron wave functions in correlated 
basis functions theory. We study the differences between the S- and D-wave 
components of the deuteron and those of the deuteron-like pair in the nuclear 
medium. The total number of deuteron type pairs is computed and com- 
pared with the predictions of Levinger's quasideuteron model. The resulting 
Levinger's factor in nuclear matter at equilibrium density is 11.63. We use 
the local density approximation to estimate the Levinger's factor for heavy 
nuclei, obtaining results which are consistent with the available experimental 

data from photoreactions. 
PACS number (s): 

Typeset using REVT^ 
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I. INTRODUCTION 



The suggestion that the nuclear response may be interpreted as the response of a collec- 
tion of neutron-proton (np) pairs carrying the quantum numbers of the deuteron was first put 
forward in the fifties by Levinger and Gottfried 0, to explain nuclear photoabsorption 
data. 

The basic idea underlying the Levinger's quasideuteron (QD) model is that the nuclear 
photoabsorption cross section aA^E^), above the giant dipole resonance and below the pion 
threshold, is proportional to that corresponding to the break-up of a deuteron embedded in 
hadronic matter, and denoted hereafter as aqD^E^) 

aA{E^) = Vd (Tqd{E^) . (1) 

The proportionality constant Vd has to be interpreted as the fraction of the A{A — l)/2 
nucleon-nucleon pairs, which are of QD type, and it is given by 

"Z(A-Z) 
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where A and Z denote the nuclear mass and charge and L is the so called Levinger's factor. 
Vd can be directly calculated from the ground state wave function of the nucleus with mass 
A. Since the deuteron is a bound state, Vd scales with the number of particles A. 

From Vd) the probability of finding a deuteron-like nucleon pair in a complex nucleus 
can be easily extracted. Such probability can be obtained by normalizing the number of 
QD pairs Vd to the total number of pairs, and, therefore, it is inversely proportional to the 
number of particles. The probability is zero in infinite nuclear matter, unless the nuclear 
matter wave function contains a long range order, providing a condensation of QD pairs. 

According to the Levinger's model |^ aqD^E^) is taken as the deuteron cross section 
times a damping function, of exponential form, accounting for the Pauli blocking of the final 
states available to the nucleon ejected from the QD: 

aQD{E,) = a,{E,) e'^^/^-) . (3) 



Subsequently, Laget proposed to associate agoiE^) with the transition amphtudes 
of virtual (vr + p)-meson exchanges between the two nucleon of the QD pair, leading to a 
cross section denoted af^'^^{E^). 

Both models fit reasonably well the existing photoreaction data in heavy nuclei, but 
the resulting factors, Lievi^) and Liageti^), have different phenomenological values, with 
LLaget{A) being about 20% larger than Lls^j^A). 

A generalization of the QD model was proposed by Frankfurt and Strikman f^, to ex- 
plain the production of fast backward protons in semi-inclusive processes off nuclear targets. 
According to the model of Refs. generally referred to as few nucleon correlation model, 
the structure of the nuclear wave function at short internucleon distances is dominated by 
strongly correlated multinucleon clusters. A quantitative understanding of the above reac- 
tion processes requires a microscopic calculation of the quasideuteron distribution P£)(k£)) 
in the nucleus, as a function of its momentum k/?. Moreover, the integral of PdO<-d) over 
ki), being proportional to Vd, provides an unbiased calculation of the Levinger's factor L. 

More recently, the occurrence and spacial structure of deuteron-like configurations in 
light nuclei has been studied using the Green's Function Monte Carlo (GFMC) method [Q. 
It is interesting to extend such analysis to heavier nuclei and to nuclear matter. 

Systematic quantitative investigations of nucleon-nucleon (NN) correlations in nuclear 
matter have been carried out within microscopic many-body theories (for a recent review see 
Ref. 0). In particular. Correlated Basis Function (CBF) theory has been applied to obtain 
the nuclear matter momentum distribution and spectral functions p|- |T^] from realistic 
hamiltonians. In this paper we use the same many body framework to carry out an ab initio 
calculation of the momentum distribution PD(kz)) of QD pairs in infinite nuclear matter, as 
well as of the associated total number of QD pairs per particle Vd/A. 

The definition of the QD total momentum distribution in terms of the overlap between 
the nuclear matter and the deuteron ground state wave functions is given in Sec. II, where 
the many-body formalism employed in the calculations is also briefly outlined. In Sec. Ill 
the results of numerical calculations, including both the QD momentum distribution and 
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Vd in nuclear matter at the empirical saturation density, p = 0.16 fm~^, are discussed and 
compared to the empirical estimates of the Levinger's factor. Finally, the summary and 
conclusions are given in Sec. IV. 



II. FORMALISM 



The distribution of QD pairs with total momentum k^) in nuclear matter is defined as 
(sum over repeated greek indeces is implicit hereafter) 

1 



PD(kz,) = — — -EEKr(ki^ 



(4) 



where J^^ = 1 is the spin of the deuteron, and 



M^^^iko) = I dR£n£r2^%,,ir,,r2,R)ni^i,r2)^niR) 



(5) 



with R = (r3, . . . , r^). In the above equation, "^nm and denote the normalized nuclear 
matter ground state wave function and the wave function of the (A — 2)-nucleon system in 
the state n, respectively. The configuration space deuteron wave function (DWF) can be 
cast in the form 



(6) 



where Q is the normalization volume, Rjj = (r^ + rj)/2, = r j — r^, |00) is the spin-isospin 
singlet two-nucleon state and the relative motion of the pair is described by 



rD{^3) 



(7) 



In Eq.(^, Uoir) and Woir) are the i = and i = 2 components of the deuteron wave 
function, normalized according to 



r dr 



(8) 



erf (a = 1, 2, 3) denote the Pauli matrices and the tensor operator is given by 
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T"^(f.,) = Sfr^.rg - 5"^ . (9) 

In CBF theory I^^nm) is usually written, in coordinate space, in the form (R= (ri, . . . , r^) 
specifies the nucleon positions) 

^nm{R)=S[1[F{zj)]MR) , (10) 

i<j 

where S is the symmetrization operator and $o is the Slater determinant describing a non- 
interacting Fermi gas of nucleons carrying momenta k with |k| < kp = (67r^p/z/)^/^, u being 
the degeneracy of the momentum states (in symmetric nuclear matter z/ = 4). The oper- 
ator F{ij), accounting for the correlation structure induced by the nucleon nucleon (NN) 



interaction, has been chosen of the form 13 



= fc{rij) + fa{rtj){(Ti ■ (Tj) + fr{rij){T, ■ Tj) + faT{rij){(Ti ■ CTj){T, ■ Tj) 

+ ft{r^j)Tafs{rij)a°'a^ + ftT{ry)Taf3{rij)aia^{Ti ■ Tj) , (11) 

where fc{r), /o-(r), /r(r), far^r), ft{r) and ftrir) are correlation functions whose radial shapes 
are determined minimizing the expectation value of the hamiltonian in the ground state 
described by Eq.(|TO|) ||T^. As r — > oo, /c(r) — 1, while all other correlation functions go to 
zero. 

Summation over k/j of Pc>(k/5) yields VdH'^Jd + 1) in nuclear matter. This number, 
which corresponds to an extensive quantity and therefore is propotional to A, leads to a 
direct evaluation of the Levinger's factor L, to be compared with the value resulting from 
the phenomenological analyses of the available experimental data on photoreactions 

1,111 • 



The quantity defined by Eqs.(^ and (^ is related to the fully linked part of the two- 
nucleon density matrix, p*^^^ (ri, r2, ri/, r2'). This part is the only one providing extra informa- 
tion on the N-N correlations with respect to that carried by the one-body density matrix, or, 
equivalently, by the nucleon momentum distribution [0. Using standard cluster expansion 
techniques [0, Pd{}^d) can be written as a series of terms involving an increasing number 
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of particles. We have calculated the cluster contributions associated with the diagrammatic 
structure shown in fig.|l|, and its exchange counterpart, where the deuteron wave function 
\1/d(1,2) is multiplied by the correlation operator F{1,2). This corresponds to a dressed 
leading order approximation, whose validity has been checked in previous CBF calculations 
of the response function and of the spectral function of nuclear matter, and whose expression 
reads 

Poi^^D) = \^j d^rn'd'rud^rr2' e^kz, (r,„+r,„)/2^ (^^^,) ^ ^(ris, rr2') , (12) 

with 

S(ri2, rr2') = ^Tr [Ft(l'2')^S^(l'2')nooV^S(12)i^(12) (1 - P.P.)] . (13) 

In the above equation, IIoo is the operator projecting onto the S = 0, T = two-nucleon 
state: 

noo - ^ ' • '-<^;-"^' . (14) 

while the spin- and isospin-exchange operators, P^r and P., are given by 

P^ = l±i^ll^^ P. = l±ilil^. (15) 



The function n{r) is the correlated one-body density matrix [§, normalized as n(r = 
0) = 1, and trivially related to the nucleon momentum distribution, n{k), through 

n(r) = I d'M^Mk) . (16) 

Evaluation of the trace appearing in Eq.(|13D leads to the simple result 

E(r, r') = ^ [U{r)U{r') + W{r)W{r')Q{v, r')] , (17) 

where 

g(r,f') = ^[3(r.rO'-l] , (18) 
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and the functions U{r) and W{r) are defined as 

U{r) = U£){r) — Au{r) , 



(19) 



W{r) = WEiij) — Aw{r) . 



(20) 



The exphcit expression of the functions Au{r) and Aw{r), yielding the deviation of U{r) 
and W{r) from the bare components of the DWF, are 



Au{r) = uoir) [K{r) - /,(r) + 3/,(r) + 3/,,(r)] 
- y/8wn{r) [/,(r) - 3fUr)] , 



(21) 



and 



Awir) = WDir) [h,ir) - ^(r) + 3/,(r) + 3^,(r)] 
'uo{r) - ^) [Mr) - 3fUr)] , 



(22) 



with hc{r) = 1 — fdr). Note that in absence of correlations, i.e. setting fdr) =1 and all 
other correlation functions identically equal to zero, U{r) and W{r) reduce to U£){r) and 
woir), respectively. 

Using the functions defined in Eqs.(|l^) and (0), the wave function describing the motion 
of the QD pair in nuclear matter can be written in the same form as the DWF (see Eqs.(|^), 

(I) and my. 

^^^J^T°^ff)afl |00) . 



V2 



(23) 



Using the Fourier tranforms of U{r) and W{r), defined as 

^ /-oo 



U{k) 



TT JO 



r drjo{kr)U{r) 



(24) 



and 



W{k) 



TT JO 



r"^ dr i2{kr)W {r) 



(25) 



jo{kr) and j2{kr) being spherical Bessel functions, and the nucleon momentum distribution 
in nuclear matter, n{k), Eq.([I^) can be rewritten in the form 



where 



P(kz5,k) 



n 



2 



n 



2 



and 



\^QD{k)\' 



(26) 



(27) 



i^QDm' =^[u\k) + w\k) 



(28) 



The above equations have been used to carry out the numerical calculations. 

It has to be noticed that the contributions arising from the non commuting structure of 
the correlations reaching the four external vertices, 1, 2, 1' and 2', of the diagrammatical 
structure of fig. |l|, are not exactly accounted for, but only according to the dressed leading 
order approximation. 



III. RESULTS 

Fig. H shows the behavior of U{r) and W{r) evaluated using a many body hamiltonian 
including the Urbana Vi4 NN potential and supplemented by the TNI model of many-body 



forces [T^. For comparison, we also show the components of the Urbana fi4 DWF and 
the functions Am and Aw defined in Eqs. (|21|) and (0), respectively. It appears that the 
main differences between deuteron and QD occur at r < 2 fm. At small relative distance 
(r < 1 fm), the effect of the nuclear medium leads to an appreciable suppression of U{r) 
with respect to uoir), whereas Wn^r) turns out to be substantially enhanced, compared to 

The momentum space behavior of \U{k)\, \W{k)\, \uD{k)\, \wD{k)\, \Au{k)\ and \Aw{k)\ 
is displayed in fig. ^ The main effect of the nuclear medium appears to be a shift of the 
second mimimum of both |f/(A;)| and |iy(fc)| towards lower values of k. 
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Eqs.(pTD and (p^) show that the nuclear medium modifications to the DWF are driven 
by the functions Ht{r) = ft{r) - 3ftr{r) and AHc{r) = -fa{r) + 3/^(r) + 3f„r{r), resuhing 
from the combination of different components of the NN correlation operator. The radial 
dependence of Ht{r) and AHc{r), illustrated in fig. ^, shows that the effect of scalar and 
spin-isospin correlations, described by AHc{r), dominates at very short relative distance, 
whereas Ht{r), accounting for tensor correlations, has a significantly longer range. 

The distribution of deuteron pairs with total momentum ho, PdO<-d), resulting from 
our approach is displayed in fig. ^ as a solid line. Within the Fermi gas model, PdO<-d) = 
at |k/)| > 2kF, implying that the high momentum tail of Poi^n) is entirely due to NN 
correlations. The distribution of deuterons in a Fermi gas is represented by a dashed line 
in the figure. The comparison between the two curves clearly shows that the correlations 
deplete the distribution with respect to the Fermi gas at |k£)| < 2kp. The depletion is 
mostly due to the non central, tensor correlations. 

Similarly, one can define the relative momentum distribution of the nucleons belonging 
to a QD pair in nuclear matter 

p^^'{k) = m\^QD{k)\\ (29) 

where 




For example, in the Fermi gas model n{k) = 6{kp — k), and (pik) takes the simple form 
m = {2Tif2p + ]^x'^ e{l - x) , (31) 

with X = k/kp. 

Fig. ^ shows the relative momentum distribution of a QD pair in nuclear matter, as 
well as the functions 1^^1315(^)1^ and (pik) defined by Eqs.(p8|) and (pOf), respectively. For 
comparison the relative momentum distribution of a deuteron in free space is also displayed. 
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The total number of pairs of the QD type in nuclear matter, Vd^ can be obtained by 
momentum integration of either P^'^'(k) or PD(k£)) times the spin multiplicity, 27^1 + 1 = 3, 
of the deuteron: 

- ^ ^ ^'^^ Pd{\^d) = - / 7^ Pn\^) ■ (32) 



A p J (27r)3 p J (27r)= 

The calculation carried out using the correlated model of nuclear matter and Eq.(p6D yields 
Vd/A = 2.895, to be compared with the Fermi gas model result of 3.406. 

In order to compare the calculated Vd to the number of QD pairs extracted from the 
analysis of photonuclear data we have to make a connection with the Levinger's formula 
given in Eqs.(|]) and (j^). The relation is given by 

P„ . L . (33) 

and, for symmetrical matter [Z = A/2), one has 

L{A) = 4 ^ . (34) 

The nuclear matter value resulting from our calculation gives L(oo) = 11.63. This value 
should be compared with that given by the phenomenological formula 

= 13.82^ , (35) 

reported in Ref. [jl4|, providing L^evioo) = 9.26. Notice that, for a deuteron in a Fermi 
gas, Lfg{oo) = 13.6. Surface contributions to L{A) can be estimated by exploiting the 
calculation of the enhancement factor JC in the electric dipole sum rule for finite nuclei of 
Ref. [^, performed within the CBF theory and Local Density Approximation (LDA). The 



enhancement factor is related to experimental data on photoreactions through the equation: 



1 ~l~ ^exp 



/rriTrC' 
aA{E^)dE^ , (36) 



where cxo = 60 [Z{A — Z)/A] MeV mb and 171^^0^ is the vr-meson production treshold. 

Therefore, the Levinger's factor can be related to /C in the mass number range where 
the coefficient D in Eq.(^ is fairly A-independent, namely for sufficiently large values of A. 
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By adding the surface contributions, as extracted from Ref . , to the nuclear matter bulk 
result, we get: 

L{A) = 11.63 - 9.76 A'^^^ . (37) 
Fig. 1^ shows our results for L{A) compared with L^evi^) and Liageti^), as extracted 



l^ , p!5|l from the available experimental data on photoreactions. The computed Levinger's 
factors are almost A-independent for heavy nuclei {A > 100), and result to be ~ 25% 
larger than LieviA) and ~ 15% smaller than LiagetiA), and therefore they are consistent 
with the experimental data. This differs from what happens for the enhancement factor /C, 
where essentially the same theory as the one used in this paper leads to a value which is 
~ 60% larger than the experimental one. Therefore, this disagreement between theory and 
experiment has to be mainly traced back to the sizeable tail contributions to the electric 
dipole sum rule, absent in the definition of Eg. (|36D . 



IV. CONCLUSIONS 

Correlated Basis Function theory of the two-body density matrix has been applied to 
compute the distribution P£)(k£)) of neutron-proton pairs characterized by the deuteron 
wave function and having total momentum k^). 

It has been found that this distribution in nuclear matter is mostly concentrated at 
< |kz)| < 2kF- Besides being responsible for the appearance of the tail of Pd(}<-d) at 
|k£)| > 2kF, NN correlations produce an appreciable effect at low momenta. The inclusion 
of correlations associated with the tensor component of the one pion exchange interaction 
leads to a ~ 15% decrease of PD^kn) at |ki5| < kp. In general, inclusion of correlations 
reduces the prediction of the Fermi gas model in this region . 

Summation of Poikn) over k/) provides the total number Vd of QD pairs, and, conse- 
quently, allows for an ab initio calculation of the Levinger's factor, L{A). The CBF results 
for symmetrical nuclear matter, L(oo) = 11.63, is about 20% larger than L^evi'^) = 9.26 
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given in the literature. In the case of heavy nuclei L^eviA) and LiagetiA) bracket our 
CBF results, which are therefore consistent with the available photoreaction data within the 
quasideuteron model phenomenology. 

It should be noticed that the theoretical estimate of L{A) in the range 150 < A < 250 
is fairly constant, its increase with A being of ~ 3%. The A~^^^ surface behavior leads to 
a very slow increase of L{A) with A, and at A ~ 200 we are still quite far away from the 
asymptotic region. 

In addition, the analysis described in this paper shows that when a deuteron is embedded 
in nuclear matter at equilibrium density, its wave function gets appreciably modified by the 
surrounding medium. While in the case of the S-wave component the difference is mostly 
visible at small relative distance (r < 1 fm), the D-wave component of the QD appears to 
be significantly quenched, with respect to the deuteron wc^r), over the range < r < 2 
fm. It has to be pointed out, however, that the radius of the QD configuration is very close 
to the deuteron radius, the difference being ~ 2 %. This result is in agreement with the 
conclusions of a recent study of deuteron-like configurations in light nuclei |^ . The authors 
of Ref. 1^ find that the density distributions of np pairs carrying the deuteron quantum 
numbers in '^He, ^He, ^Li, ^Li and ^^O exhibit size and structure similar to those observed 
in the deuteron. 

The relative momentum distribution of a QD pair, PJ^\k), extends into the region |k| > 
kp, where it appears to be strongly suppressed with respect to the corresponding deuteron 
momentum distribution |\l/£)(A;)p, although |\1/q£)(A;)P is larger than |\I'£)(A;)p at high k. It 
has to be pointed out that the behavior of P^''{k) aX k > kp is entirely dictated by the high 
momentum tail of the nuclear matter momentum distribution, produced by strong short 
range NN correlations. Within the Fermi gas model n{k > kp) = 0, and P^^^k > kp) 
vanishes identically. 

Higher order cluster terms, neglected in this paper and arising from the inclusion of addi- 
tional bonds in the diagrammatical structure of fig. |I|, are not expected to change the main 
conclusions of the present paper, neither regarding the behavior of the deuteron distribution 
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in nuclear matter, nor as far as the disussion on the Levinger's factor is concerned. 

In view of the relevance that PD^ko) and |\['q£i(A;)P may assume in the study of those 
lepton-nucleus reactions where the ejected hadron is in kinematical regions forbidden to 
lepton-nucleon processes, the calculations presented in this paper need to be extended i) 
by introducing higher order cluster terms in the expansion of the two-body density matrix, 
and ii) by explicitely considering finite nuclei wave functions. Work in these directions is in 
progress. 
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FIGURES 

FIG. 1. Diagram showing the cluster contribution to Pd{^d) of eqs.(^) and (|5|) considered in 
this paper. The oriented solid lines represent the correlated one-body density matrix, whereas the 
wiggly lines correspond to the dressed deuteron-like np pairs. 

FIG. 2. Upper panel: the solid line shows the radial dependence of U{r), defined by Eq.(|l9|), 
while the dashed and dot-dash lines correspond to UD{r) and An(r), respectively. Lower panel: 
same as the upper panel, but for the i=2 components of the QD and deuteron wave functions. All 
wave functions are given in units of (GeV/c)^/^. 

FIG. 3. Same as in fig. ^ in momentum space. All wave functions are given in units of 
(GeV/c)-3/2. 

FIG. 4. Radial dependence of the functions AHc{r) and Ht{r), entering the definitions of 
Au{r) and Aw{r) (see Eqs.(|TD and (|2|). 

FIG. 5. Momentum distribution of QD pairs in nuclear matter at equilibrium density as a 
function of the total momentum |k£)| (see. Eqs.(^ and (|5|)). Solid line: correlated model; dashed 
line: deuterons in a Fermi gas model. The insert shows a blow up of the region \]i.£)\/2kF < 1, 
plotted in linear scale. 

FIG. 6. The solid line shows the relative momentum distribution of a QD pair in nuclear 
matter at equilibrium density (Eq.(pg|)). The dashed and dot-dash lines correspond to (l){k) (in 
units of (GeV/c)3) and |^'qd(A;)|2 (in units of (GeV/c)-^), defined by Eqs.(|oD and (H). The 
diamonds show the squared momentum space wave function of a free deuteron. 

FIG. 7. The CBF Levinger's factor L{A) of heavy nuclei (solid line) and nuclear matter 
(indicated by the arrow). The LDA approximation of ref. has been used for heavy nuclei. The 
phenomonological values of L^eviA) corresponding to photoreaction data of Lepretre et al. ||l^ 
(squares) and Ahrens et al. (crosses and diamonds) are taken from ref. \14\. The empirical 
values of L Lev (A) represented by circles are from ref. ||2l| . 
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